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ABSTRACT 

Aims. Asymptotic giant branch (AGB) stars are among the largest distributors of dust into the interstellar medium, and it is therefore 
important to understand the dust formation process and sequence in their strongly pulsating extended atmosphere. By monitoring the 
AGB star W Hya interferometrically over a few pulsations cycles, the upper atmospheric layers can be studied to obtain information 
on their chemical gas and dust composition and their intracycle and cycle-to-cycle behavior. 

Methods. Mid-infrared (8-13 jjm) interferometric data of W Hya were obtained with MIDI/VLTI between April 2007 and 
September 2009, covering nearly three pulsation cycles. The spectrally dispersed visibility data of all 75 observations were ana- 
lyzed by fitting a circular fully limb-darkened disk (FDD) model to all data and individual pulsation phases. Asymmetries were 
studied with an elliptical FDD. 

Results. Modeling results in an apparent angular FDD diameter of W Hya of about (80 ± 1.2) mas (7.8 AU) between 8 and 
10 yum, which corresponds to an about 1.9 times larger diameter than the photospheric one. The diameter gradually increases up 
to (105 ± 1.2) mas (10.3 AU) at 12 /im. In contrast, the FDD relative flux fraction decreases from (0.85 ± 0.02) to (0.77 ± 0.02), 
reflecting the increased flux contribution from a fully resolved surrounding silicate dust shell. The asymmetric character of the ex- 
tended structure could be confirmed. An elliptical FDD yields a position angle of (11 ± 20)° and an axis ratio of (0.87 ± 0.07). A 
weak pulsation dependency is revealed with a diameter increase of (5.4 ± 1.8) mas between visual minimum and maximum, while 
detected cycle-to-cycle variations are smaller. 

Conclusions. W Hya's diameter shows a behavior that is very similar to the Mira stars RR Sco and S Ori and can be described by an 
analogous model. The constant diameter part results from a partially resolved stellar disk, including a close molecular layer of H2O, 
while the increase beyond 10 yum can most likely be attributed to the contribution of a spatially resolved nearby AI2O3 dust shell. 
Probably due to the low mass-loss rate, close Fe-free silicate dust could not be detected. The results suggest that the formation of 
amorphous AI2O3 occurs mainly at visual minimum. A possible close AI2O3 dust shell has now been revealed in a few objects calling 
for self-consistent dynamic atmospheric models including dust formation close to the star. The asymmetry might be explained by an 
enhanced dust concentration along an N-S axis. 

Key words, stars: AGB and post-AGB - stars: individual: W Hya - stars: circumstellar matter - stars: diameter ~ infrared: stars - 
techniques: interferometric 

1. Introduction The M-type (O-rich) star W Hya is interferometrically ob- 

served in the thermal N-band to resolve the upper atmospheric 
Asymptotic giant branch (AGB) stars have high luminosities, a layers over a few pulsations cycles. This gives the possibility 
high mass-loss rate, and a relatively cool atmospheric temper- of obtaining information on its chemical gas and dust composi- 
ature, and they represent a late stage of stellar evolution. Most tion and its behavior throughout and between pulsation cycles, 
of the AGB stai's pulsate (as Miras, semi-regular valuables or This can lead to a better understanding of the complex dynam- 
irregular variables) and have typical main-sequence progenitor ical atmospheric processes attributed to AGB stars. Since AGB 
^ masses of 0.8 to 10 Mq. A large fraction of the mass is con- stai's are one of the most important distributors of dust into the 
centrated in a tiny degenerate carbon-oxygen core, surrounded interstellar- medium, besides red super giant (RSG) stars and su- 
■ ■ ■ by a helium- and hydrogen-burning shell and an extended stellar pernovae (SN), it is especially important to understand the dust 
envelope with a very extended convective zone. Photospheric di- formation process and sequence in their strongly pulsating ex- 
ameters are typically a few AU. The star is embedded in a dusty tended atmosphere, 
circumstellar envelope (CSE) (cf. e.g. review book by Habing & 

Olofsson.2004il. t jj- ■ • . ^, 

In addition, asymmetries (e.g. oblateness due to rotation or 

a companion, or non-symmetric brightness distributions due to 



* Based on observations made with the Very Large Telescope stellar spots or large convection zones), which are not uncom- 

Interferometer (VLTI) at the Paranal Observatory under program IDs for AGB stars, can be investigated. Past asymmetry deter- 

079.D-0140, 080.D-0005, 081.D-0198, 082.D-0641 and 083.D-0294. minations of W Hya were not very conclusive, in the sense that 

** Color versions of the figures are available in electronic form via position angles over a wide range were repoited, ranging from 

http://www.aanda.org about 70° (Szymczak et al. 1998) in the radio to about 143° 

Fellow of the International Max Planck Research School (IMPRS), (Lattanzi et al. 1 19971 1 in the visual, while no departures from 



*** 



e-mail: rgeisler(aisw. uni -Heidelberg . de symmetry could be detected in the near-infrared (IR) within the 
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measurement uncertainties (Ireland et al. I2004al and Monnier et 
al. 120041) . 

W Hya is one of the best observed AGB stars in the southern 
hemisphere. In particular, many interferometric diameter mea- 
surements at visual and near-IR wavelengths have been car- 
ried out (e.g. Ireland et al. 2004a and Woodruff et al. 120091 
cf. Sect. 13.3b . Because of the very extended atmosphere, it is 
expected that different wavelengths probe different atmospheric 
layers in AGB stars (Baschek et al. [TWT] and Scholz et al. l200TT l. 
The visual and near-IR diameters of the star W Hya range from 
about 30 up to 70 mas, correlated to the absorption and emission 
bands of the most abundant and radiatively important molecular 
species, such as H2 O, OH , CO, TiO, SiO, CO2 and SO2, besides 
H2 (Hofmann et al. [19981 and Jacob et al. 2000). Therefore it is 
difficult to measure a continuum diameter, while these opaque 
molecular layers can easily have radii twice the continuum ra- 
dius (e.g. Mennesson et al. 120021 Tej et al. 120031 and Ohnaka 
12004 !) with temperatures of 1000 to 2000 K, i.e. below the tem- 
perature of the continuum-forming surface of the star, but higher 
than the temperature of the surrounding circumstellar dust enve- 
lope. 

Dust grains with high sublimation temperatures are already 
present in the outer atmosphere and are the seed particles for 
further dust growth (e.g. Lorenz-Martins & Pompeia 2000 1 and 
Verhoelst et al. 120091 and references therein). Radiative pres- 
sure leads then at larger radii to an acceleration of the dust and 
the gas (through frictional coupling), resulting in the character- 
istic high mass-loss rates of AGB stars (e.g. Woitke 120061 and 
H6fner^2008 ). Several attempts have been made to include dif- 
ferent dust species, such as aluminum oxide (AI2O3), silicates 
(e.g. Mg2Si04), spinel (MgAl204), and olivine (MgFeSi04), to 
model the spectral energy distribution and to explain the nearly 
two times larger mid-IR diameters by adding a dust shell to a 
dust-free atmospheric model (cf. e.g. for Mira stars Ohnaka et 
al. 12005 and Wittkowski et al.l2007l and for RSG stars Perrin et 
al. l2007i and Verhoelst et al. 120091 ). The radiative importance of 
dust in an extended molecular shell depends on the local tem- 
perature, the time dependent growth rate, and the wavelength- 
dependent optical properties. The former can be very different 
for different pulsation phases and cycles. 

Dynamical atmospheric and wind models for AGB stars 
were developed by Hofmann et al. ( 119981 ). Woitke et al. (119991 
I2OO6I ). Scholz et al. (|200TT ). Hof ner et al. (12003). Ire land & 
Scholz '20061 Ho fner & Andersen (12007^ . Ire land e t al. (l2008l ). 
Nowotny et al. (120101 ) and Lebzelter et al. (120101 ). The exact 
chemical composition, pulsation phase, and pulsation cycle be- 
havior of AGB stars, e.g. the location of molecular and dust lay- 
ers, is still very uncertain. Observations presented in this work 
can help to improve this situation and provide the opportunity to 
compare with theoretical model predictions. 

W Hya is a large-amplitude, semi-regular variable (SRa, but 
sometimes also classified as Mira) with a pulsation period of 
about one year, and it is located in the P-L-diagram on se- 
quence C (fundamental mode pulsator, Lebzelter et al. 120051 
sequence 1 in Riebel et al. '2010). The visual magnitude varies 
strongly, while the amplitude in the N-band is rather small 
(cf. Sect. 12.31 ). The radial velocity amplitude, derived from CO 
Av = 3 lines, is about 15 km s ' (Hinkle et al. 119971 Lebzelter et 
al. l2005l ). Distance estimates to W Hya range from 78 pc (Knapp 
et al. 120031 revised Hipparcos value) to 115 pc (Perryman et 
al. 119971 Hipparcos value). Throughout this paper an interme- 
diate value of 98;'^jg pc from Vlemmings et al. (120031 ) will be 
assumed (Very Long Baseline Interferometry maser measure- 



ment). W Hya's luminosity is about 5400 (Justtanont et 
al. l2505l ). 

After a description of the observational method and data re- 
duction in Sect. |2] the model is explained in Sect. [3] where the 
results will be interpreted by the presence of dust close to a 
molecular layer. Sect. U investigates the pulsation dependence 
of the apparent diameter, as well as departures from symmetry. 
A summary is given in Sect.|5] 

2. Observations and data reduction 

2.1. interferometric observations witti l\/iiDi/VLTi 

The data presented here were obtained with the mid-IR (8 - 
13 yum) interferometer MIDI (Leinert et al. l2003l 12004^ at 
the Very Large Telescope Interferometer (VLTI) in service 
mode. W Hya is one of five stars monitored from P79 to P83 
(April 2007 to September 2009) under the program IDs 079. D- 
0140, 080.D-0005, 081.D-0198, 082.D-0641, and 083.D-0294. 
An observation log is given in Table[T] showing all 83 individual 
observations. Six different baseline configurations of two auxil- 
iary telescopes (ATs) were used. This results in projected base- 
lines from 13 to 71 meter and a wide range of position angles 
(PA, §; east of north). The uv-coverage of all used data can be 
seen in the right hand panel of Fig. [1] 

Before or after each target observation the calibrator 2 Cen 
(HD 120323) was observed with the same setup to calibrate the 
visibility measurements. The calibrator star has a diameter of 
(13.25 ± 0.06) mas (model diameter from Verhoelst 200^, a 
specti-al type of M4.5 III, a 12 yum flux of (256 + 26) Jy (IRASB, 
and an angular separation to the target of about 6°. Observations 
were executed in SCI -PHOT mode, where the photometric and 
the interferometric spectra are recorded simultaneously. This has 
the advantage that the photometry and the fringe signal are ob- 
served under the same atmospheric conditions. Either the prism, 
with a spectral resolution of R - A/ AA - 30, or the grism, with 
a spectral resolution of 230, were used to obtain spectrally dis- 
persed fringes. 

2.2. MiDi Sci-Piiot data reduction 

The standard MIA+EW^ (version 1.6) data reduction package 
with additional routines for processing SCI -PHOT data (Walter 
Jaffe, private communication) was used. In any observation in 
SCI-PHOT mode, all four read-out windows (channels) are il- 
luminated simultaneously. The two central channels record the 
interferometric signals, while the two outer channels are dedi- 
cated to measuring photometric fluxes (cf. Fig. |2). To measure 
these signals, the two telescope beams (one beam from each 
telescope) are separated by a beamsplitter. One of the separated 
beams of each telescope is sent directly to the photometric de- 
tector window, while the other two remaining beams (one from 
each telescope) are combined in the beam combiner to obtain 
two interferometric signals with opposite phases. 

To be able to calculate the instrumental visibility by dividing 
the interferometric fringe signal (correlated flux) by the photo- 
metric signal (uncorrelated flux), the exact light split ratios be- 
tween these channels have to be known as functions of wave- 
length. This requires recording additional photometry. By clos- 

' http://www.ster.k;uleuven.ac.be/~tijl/MIDLcalibration/mcc.txt; see 
also ESO CalVin database: http://www.eso.org/observ ing/etc/ 
- http://irsa.ipac.caltech.edu/Missions/iras.html 
' http://www.strw.leidenuniv.nl/~koehler/MIDI 



Table 1. Observation log. 
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Date 


Al 


jjisp 


a (m) 


fA ( ) 


Phase 




Date 


/\i 


r^i f t-»b 
uisp 


a (m) 


m (, ) 


Phase 




2007-04-12 


A 


grism 


14.54 


51.81 


0.30 


used 


2008-04-03 


B 


prism 


29.19 


82.75 


0.22 


used 


2007-04-13 


B 


grism 


29.10 


51.92 


0.30 


used 


2008-04-28 


D 


prism 


59.43 


54.70 


0.29 


used 


2007-04-17 


B 


grism 


22.34 


94.12 


0.31 


used 


2008-04-28 


E 


prism 


70.58 


130.03 


0.29 


n. u. 


2007-04-22 


D 


grism 


59.53 


54.92 


0.33 


used 


2008-04-28 


E 


prism 


71.32 


10.41 


0.29 


n. u. 


2007-04-22 


D 


grism 


59.53 


54.92 


0.33 


used 


2008-04-28 


D 


prism 


45.23 


93.72 


0.29 


used 


2007-04-22 


D 


grism 


42.06 


96.19 


0.33 


used 


2008-05-25 


D 


prism 


61.21 


58.90 


0.36 


used 


2007-04-22 


D 


grism 


38.56 


99.06 


0.33 


used 


2008-05-30 


D 


prism 


30.02 


107.46 


0.37 


used 


2007-04-24 


E 


grism 


67.30 


125.04 


0.33 


n. u. 


2008-07-03 


A 


prism 


13.67 


86.31 


0.46 


used 


2007-04-25 


E 


grism 


64.92 


123.16 


0.34 


n. u. 


2008-07-03 


C 


prism 


32.03 


95.67 


0.46 


used 


2007-04-25 


F 


grism 


71.47 


2.19 


0.34 


used 


2008-07-06 


D 


prism 


36.87 


100.52 


0.46 


used 


2007-06-18 


A 


grism 


13.25 


87.72 


0.47 


used 


2009-01-16 


F 


prism 


71.23 


-12.23 


0.96 


used 


2007-06-20 


C 


grism 


46.98 


62.82 


0.48 


used 


2009-01-16 


D 


prism 


60.22 


56.51 


0.96 


used 


2007-06-20 


C 


grism 


43.73 


82.84 


0.48 


used 


2009-01-17 


E 


prism 


65.19 


123.34 


0.97 


n. u. 


2007-07-02 


c 


grism 


48.00 


70.64 


0.51 


used 


2009-01-20 


C 


prism 


46.78 


61.99 


0.97 


used 


2007-07-02 


c 


grism 


40.49 


86.86 


0.51 


n. u. 


2009-01-21 


B 


prism 


25.33 


32.67 


0.98 


used 


2007-07-04 


B 


grism 


31.55 


75.51 


0.52 


used 


2009-01-21 


B 


prism 


28.19 


47.79 


0.98 


used 


2007-07-04 


B 


grism 


24.39 


90.99 


0.52 


used 


2009-01-21 


B 


prism 


29.88 


55.45 


0.98 


used 


2008-01-10 


B 


prism 


30.08 


56.39 


0.01 


used 


2009-01-22 


C 


prism 


43.10 


50.26 


0.98 


used 


2008-02-20 


D 


grism 


59.95 


55.89 


0.11 


used 


2009-01-22 


C 


prism 


44.68 


55.03 


0.98 


used 


2008-02-20 


D 


prism 


62.76 


63.21 


0.11 


used 


2009-01-22 


A 


prism 


45.94 


59.00 


0.98 


used 


2008-02-21 


D 


prism 


63.89 


68.30 


0.11 


n. u. 


2009-01-25 


A 


prism 


12.34 


28.13 


0.99 


used 


2008-02-22 


B* 


prism 


31.85 


73.49 


0.12 


used 


2009-01-25 


A 


prism 


13.05 


37.34 


0.99 


used 


2008-02-22 


B* 


prism 


31.00 


78.00 


0.12 


used 


2009-01-25 


A 


prism 


14.35 


50.06 


0.99 


used 


2008-03-02 


E 


prism 


66.93 


124.69 


0.14 


used 


2009-01-25 


A 


prism 


15.78 


64.41 


0.99 


used 


2008-03-02 


F 


prism 


71.45 


4.51 


0.14 


used 


2009-01-27 


C 


prism 


41.69 


45.93 


0.99 


used 


2008-03-03 


D 


prism 


63.91 


72.10 


0.14 


used 


2009-01-27 


C 


prism 


43.52 


51.53 


0.99 


used 


2008-03-06 


D 


grism 


61.71 


78.46 


0.15 


used 


2009-02-16 


D 


prism 


63.97 


69.21 


0.04 


used 


2008-03-11 


A 


grism 


12.33 


28.02 


0.16 


used 


2009-02-16 


D 


prism 


63.24 


75.18 


0.04 


used 


2008-03-11 


A 


prism 


13.98 


46.70 


0.16 


used 


2009-03-16 


D 


prism 


62.93 


76.01 


0.12 


used 


2008-03-12 


A 


grism 


15.23 


79.71 


0.17 


used 


2009-04-20 


B 


prism 


20.28 


97.36 


0.21 


used 


2008-03-13 


B 


grism 


30.54 


79.42 


0.17 


used 


2009-04-23 


C 


prism 


45.26 


80.46 


0.21 


used 


2008-03-13 


B 


prism 


31.94 


68.34 


0.17 


used 


2009-04-23 


C 


prism 


36.25 


91.34 


0.21 


used 


2008-03-13 


B 


prism 


29.95 


81.01 


0.17 


used 


2009-04-24 


B 


prism 


27.75 


85.53 


0.22 


used 


2008-03-14 


A 


prism 


15.86 


65.67 


0.17 


n. u. 


2009-05-02 


D 


prism 


55.99 


46.91 


0.24 


used 


2008-03-14 


A 


prism 


15.80 


75.37 


0.17 


used 


2009-05-03 


E 


prism 


69.19 


127.28 


0.24 


used 


2008-03-14 


A 


prism 


15.26 


79.55 


0.17 


used 


2009-05-03 


F 


prism 


71.42 


7.00 


0.24 


used 


2008-03-25 


C 


prism 


47.86 


67.68 


0.20 


used 


2009-05-03 


D 


prism 


52.46 


38.03 


0.24 


used 


2008-03-25 


C 


prism 


47.38 


75.34 


0.20 


used 


2009-06-04 


C 


prism 


45.75 


58.37 


0.32 


used 


2008-03-25 


c 


prism 


41.85 


85.29 


0.20 


used 


2009-06-04 


C 


prism 


46.29 


60.19 


0.32 


used 


2008-04-01 


B 


prism 


28.32 


48.37 


0.22 


used 


2009-06-04 


A 


prism 


12.40 


90.40 


0.32 


used 


2008-04-02 


C 


prism 


44.91 


81.06 


0.22 


used 


2009-08-15 


A 


prism 


11.73 


92.45 


0.51 


used 


2008-04-02 


A 


prism 


12.96 


88.65 


0.22 


used 

















" AT stations: A = EO-GO, B = GO-HO, B* = AO- DO, C = EO-HO, D = DO-HO, E = DO-Gl and F = HO-Gl; dispersive element; " projected 
baseline length; position angle of the projected baseline on the sky; quality flag showing if that observation is used for the model fitting or not 
(n. u.), see Sect. l3.1l for reasons that a value had not been used 



ing the shutter of one telescope, the light from the other tele- 
scope illuminates one of the photometric channels and both in- 
terferometric channels. From these data, the light split ratio be- 
tween photometric and interferometric channels for this tele- 
scope beam can be computed. This is repeated for the other tele- 
scope beam as well. In principle, this would only be necessary 
once per night or less often, but is done for every observation to 
check for consistency. 

Accordingly, the first step in the reduction is the calculation 
of these ratios (or better the mapping behavior between interfer- 
ometric and photometric channels). After unchopping the addi- 
tional photometry data, a point spread function is fitted perpen- 
dicular to the wavelength direction as a function of wavelength 
to each of the three read-out windows to obtain the spectra. This 
step also includes the determination of the sky and instrumental 



background, remaining after the unchopping, and its removal. 
The spectra of all channels are then divided by each other to de- 
rive the respective split ratios. One of the largest eiTor sources 
here is the determination of the coiTect background, since it can 
result in too high or too low fluxes in affected channels and can 
therefore lead to errors in the final visibility and photometric 
spectra. Typical final errors are on the order of 10% and 30%, 
respectively. 

With the obtained light split ratios (mapping behavior), the 
SCI -PHOT data are reduced in the next step. After unchopping 
the SCI -PHOT data, the simultaneously recorded photometric 
spectra (in the photometric channels) are multiplied by the ap- 
propriate split ratios; i.e., each photometric spectrum is mapped 
to each interferometric channel to get an estimate of the uncor- 
related (photometric) fluxes at the location of the interferometric 
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wavelength X f^tm) u (l/arcsec) 

Fig. 1. Left: All 75 calibrated visibilities as function of wavelength obtained with MIDI, color-coded by visual light phase. Errors 
are omitted for clarity but given in Table |2] Right: UV-coverage of all 75 used interferometric observations. The visibility spectra 
are binned into 25 wavelength bins. 




Fig. 2. MIDI setup in SCI -PHOT mode (adopted from Leinert et 
al. l2003T l. 



channels. The average of the summed images of each channel 
gives the final raw photometry. To obtain the visibility, the inter- 
ferometric fringe signal (correlated flux) in each interferometric 
channel is divided by the estimated geometric mean image of 
the mapped uncorrected (photometric) spectra. Prior to the di- 
vision, the fringes are coherently averaged after determining the 
group delay, similar to the HIGH- SENS reduction (cf. e.g. Ratzka 
L2005). Both visibility spectra are then averaged to get the final 
raw visibiUty. A more detailed description can be found in Zhao- 
Geisler ( l20T0l l. 



These steps are taken separately for the target (W Hya) and 
the calibrator (2 Cen). With the known calibrator diameter the 
raw target visibility is calibrated with the calibrator's transfer 
function, assuming a wavelength-independent uniform disk di- 
ameter. The raw photometry is calibrated by assuming that the 
spectrum of the calibrator looks close enough to the Rayleigh- 
Jeans part of a blackbody and by using the known 12 jjm flux of 
the calibrator. The Rayleigh- Jeans assumption is not fully valid 
between 8 and 9 jim since 2 Cen shows a weak SiO absorption 
feature in its ISO spectrum. 

After reducing and calibrating all data, results with unphysi- 
cal visibilities (due to bad environmental conditions or failure of 
the reduction process) are rejected. Eight of the 83 observations 
were thus excluded and are flagged as "not used" in Table [1] 
Since several observations of W Hya were executed in some 
nights, more than one calibrator observation was available. For 
those nights, the calibrators nearest in time were used to calcu- 
late a median calibrated visibility with the standard deviation as 
a typical error. 

Grism data were interpolated to the prism grid. The wave- 
length range between 8 and 12 /im was binned into 25 wave- 
length bins to allow a faster computation of the model fits. 
Within each wavelength bin a mean visibility and error were 
calculated. For the fits later on, we always checked that this ap- 
proach does not mask any additional spectral features, with the 
result that in all cases the shapes of the wavelength-dependent 
parameters were not altered significantly. Measurements beyond 
12 jjm were not used because the recorded flux was too low, 
making the reduction inconsistent. 

To avoid the problem of underestimating the error and to al- 
locate an error to visibilities where only one calibrator observa- 
tion exists, all available visibility errors within the same wave- 
length bin were averaged. The resulting mean error was then as- 
signed to all visibilities within this wavelength bin. The final 75 
visibility curves are shown in the left hand panel of Fig. [T] The 
central wavelengths and assigned visibility errors of each bin are 
given in the first and second columns of Table |2] respectively. 

Assuming the same error for all measurements within a 
wavelength bin improves the model fits as well. In a fitting 
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method based on a chi-square technique (cf. Sect. 13.11 1. all visi- 
bility measurements are then represented with the same signifi- 
cance. In the case that each value has an individual uncertainty, 
values with higher errors are underrepresented, while values 
with low errors are overrepresented. Consequently, higher vis- 
ibilities, with in general higher absolute errors, are not weighted 
accordingly, and the importance of lower visibilities, at generally 
higher spatial frequencies with lower absolute errors, is overes- 
timated. Only equal weighting ensures that over the whole spa- 
tial frequency range a model fits all data well, assuming that all 
measurements are equally significant. This has the consequence 
that chi-square values are only useful in a relative sense. For this 
reason chi-square values are not given. 

2.3. Light curves 

To assign a pulsation phase to the observations, visual data from 
the American Association of Variable Star Observers (AAVSOQ) 
and the All Sky Automated Survey (ASASB Pojmanski et 
al. I2005i l are used. After binning the AAVSO data into ten-day 
bins, a simple sinusoid was fitted to the AAVSO and ASAS 
data over a period of about 10 years (2000 - 2010). The fit 
gives a period of (388 + 5) days, a mean visual magnitude of 
(7.6 + 0.1) mag, a semi-amplitude of (1.4 + 0.1) mag and a 
Julian Date of maximum brightness of (2452922 + 5) days (de- 
fined as phase 0.0). These data, including the fit, are plotted in 
Fig. Inversus time and in Fig. |4] versus visual light phase. The 
visual amplitude of the semi-regular variable W Hya is much 
lower than typical Mira variables. 

In both plots, the times of the 75 used MIDI observations 
are included as tick marks above the light curve, color-coded 
by phase in Fig. |3] and color-coded by position angle in Fig. |4] 
This allows recognizing different distributions of position angles 
throughout the pulsation cycle. It can be seen that the position 
angles are not sampled well over the pulsation period. At around 
maximum light, position angles cluster around 50°, while around 
minimum light, position angles are mostly at around 90°. Only 
for the intermediate phase is the position angle sampling more 
uniform. This will be important in Sect. |4] where asymmetries, 
intracycle variations and cycle-to-cycle variations are investi- 
gated. The regions used for this study are already shaded in both 
figures. It is also notable that the whole pulsation cycle is not 
equally sampled because the pulsation period is around one year 
(also true for the AAVSO and ASAS sampHng). 

The MIDI photometry is shown in Fig. [3] directly below the 
tick marks with the magnitude scale given on the right and as 
flux values below the magnitude plot in Fig. |4] Only the av- 
eraged fluxes between 11.5 and 12.5 /im are plotted in order 
to be comparable with observations made with other instru- 
ments (e.g. IRAS). Even though the MIDI fluxes are afflicted 
by fairly large errors, a clear phase dependence is detectable in 
the phase-folded plot. A sinusoidal fit gives a semi-amplitude of 
(510 ± 100) Jy, which is equivalent to a flux variation on the or- 
der of 20% between maximum and minimum light, with a mean 
mid-IR flux of (4.9 + 0.2) kJy. The maximum occurs after the 
visual maximum at visual phase 0.15 + 0.05. 

This phase shift is consistent with previous studies of AGB 
stars (cf. e.g. Lattanzio & Wood 120041 Smith et al. 120061 and 
references therein, and Nowotny et al. '2010(. The reasons for 
this are that most of the flux is emitted in the infrared and 
that the visual spectrum is strongly influenced by features of 

http://www.aavso.org/ 
' http://www.astrouw.edu.pl/asas/ 



the temperature-sensitive TiO molecule (Nowotny et al. 1201 Oi l. 
Similar flux variations have been reported for other AGB stars by 
Ohnaka (private communication; preliminary results for R Car, 
R Cnc and Z Pup) for MIDI AT observations, by Wittkowski et 
al. (12007) (O-rich Mira S Ori), and Ohnaka et al. (.2007j (C-rich 
Mira V Oph) for MIDI Unit Telescope (UT) observations, and 
by Karovicova (private communication; preliminary result for 
RR Aql) for both AT and UT observations. Since the photomet- 
ric accuracy of MIDI used with the ATs is very low, these values 
should be used with caution. For this reason all data are inter- 
preted with respect to the visual light curve instead of the mid-IR 
light curve. It should also be kept in mind that folding consecu- 
tive cycles might not always be appropriate since the pulsation is 
not strictly regular. The irregularity can clearly be inferred from 
the plots. 

2.4. Spectra 

The median of all calibrated photometric MIDI spectra (scaled 
with the 12 fim flux of the calibrator as described above) is 
shown in the spectral energy distribution (SED) in Fig. |5] The 
uncertainties are given by the standard deviation. Even if this 
averaging over all pulsation phases and cycles might be ques- 
tionable, it is not important here, since the spectrum is not used 
for any modeling. In addition to fully scientifically verified spec- 
tra from ISC0 and IRAS, photometric data from HIPPARCOS, 
USNO-Bl, 2MASS, and IRA^i not corrected for reddening ef- 
fects, are plotted as well. A blackbody with an effective temper- 
ature of 2300 K is overplotted as guidance. Due to the infrared 
excess, strong metallic oxide line, and molecule absorption at 
shorter wavelengths, and dust extinction, it is not expected that 
a blackbody curve fits the spectral data in an appropriate way. 

The ISO SWS spectrum is mainly dominated by absorp- 
tion bands of H2O between 2.5 - 3.0 fim (stretching mode) and 
5.0 - 8.0 fim (bending mode), and an SiO absorption band be- 
tween 8 and 9 fim (v = 1 - 0). Distinct absorption lines of CO 
at around 2.4 jjm, OH at 2.9 - 4.0 jum, CO2 at 4.25 fim and 
SO2 at 7.4 jim can be weakly seen in the spectrum as well. 
Justtanont et al. (2004) derived from temperature investigations 
that these absorptions originate in different molecular layers in 
the circumstellar shell. The OH, CO, and CO2 absorption bands 
arise mainly from a hot (about 3000 K), dense region very close 
to the stellar photosphere, where H2O is still photodissociatedby 
shocks. The H2O and a second CO2 absorption band originate 
in a layer with a temperature of 1000 K, i.e. a molecular layer 
(molsphere0) farther out. The SiO molecule absorption arises in 
the same region where the H2O molecular shell exists and where 
SiO is still not condensed in dust grains. 

Dust emission can also be identified in the ISO spectra. The 
features between 10 and 20 fim are a combination of emis- 
sion from amorphous silicates at around 10 /vm, compact AI2O3 
at around 11 jum, and MgFeO at around 19 fj.m. Justtanont 
et al. (120041) have obtained a satisfactory SED fit to all three 
emission features. They derived a lov^ total mass-loss rate of 
(3.5-8) X lO^'^Moyr"', and a lower limit to the dust mass-loss 



http://iso.esac.esa.int/ida/ 
' all from http://vizier.u-strasb g.fr/viz-bin/VizieRl 
The occurrence of molecular layers is the result of the capability 
of forming certain molecules at a specific distance from the star and 
their dilution when reaching larger distances. Such a shell exists for 
all atmospheric molecules, and only their total abundance and radiative 
properties determine whether such a molsphere can be detected or not. 
' If compared with a typical AGB star. 
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Fig. 3. The visual light curve of W Hya, covering the period of the MIDI observations. The V magnitudes are from AAVSO (10 day 
bins) and AS AS. A simple sinusoidal fit is included to determine the phases used throughout this paper. The times of the MIDI 
observations are included as tick marks above the curve and are color-coded by visual light phase. The MIDI flux at around 12 fim 
is shown directly below the tick marks with the magnitude scale given on the right. The data within the three shaded regions are 
used for cycle-to-cycle studies in Sect. 14.31 



rate for siUcates, AI2O3 and MgFeO of 1.5 x lO^'^Moyr^^ 
1.3 X 10""'Moyr"' and 2.5 X 10" "Moyr "' respectively. Hinkle 
et al. (1997) and de Beck at al. (1201 Oi l obtained similar total 
mass-loss rates of 2 x lO^^Mgyr"^ and 7.8 x lO^^Moyr"', re- 
spectively. The derived low dust mass-loss rate and pulsation 
behavior is typical of a star at the beginning of the AGB phase 
and the transformation from a semi-regular variable to a Mir a 
star (Hinkle et al. 1997). W Hya is classified in the dust emis- 
sion scheme of Sloan and Price (1998) as SE8 (classic narrow 
silicate emission type). 

If the MIDI spectrum is compared with the ISO SWS spec- 
trum, it is obvious that the silicate emission is not detected. This 
is shown in more detail in the inset of Fig. |5] In the MIDI spec- 
trum of W Hya, a weak emission of amorphous AI2O3 at around 
1 1 /vm is present, while no silicate feature at around 10 jum can 
be seerPl^ This behavior can be attributed to instrumental char- 
acteristics. ISO has a much larger field of view (FOV) compared 
to MIDI. With a small FOV of about 1 to 2 arcsecO the emis- 
sion of the extended silicate dust shell is not detected in the MIDI 
spectrum. 

However, this non-detection allows a lower limit to be de- 
rived for the inner boundary of the silicate dust shell. Assuming 
a conservative value of the FOV of 1 arcsec, the main emission of 
the silicate dust shell in W Hya originates in a region with an in- 
ner radius larger than 28 photospheric radiiP^ (>50 AU. >0.5 arc- 
sec). This is consistent with knowing that the dust envelope is 



Fluctuations between 9.3 and 10 fim are caused by difficulties in the 
reduction caused by the telluric ozone feature, and flux measurements 
beyond 12 //m are probably not calibrated well. 

' ' The exact value depends on the used AT baseline, the slit/mask po- 
sition and other instrumental specifications. 

The value of the photospheric radius, Rp],^, = 6phot/2, is defined in 
Sect. [331 



very extended, up to 40 arcsec (Hawkins 1 1990) . The higher flux 
measured with MIDI compared to ISO can be attributed to the 
circumstance that the averaged spectrum is dominated by obser- 
vations conducted around mid-IR maximum phases, while the 
ISO spectrum were obtained in a post visual minimum phase 
(phase 0.8, which is close to the mid-IR minimum; cf. Fig. |4]i. 
Together with the calibration uncertainty of the averaged MIDI 
spectrum of larger than 10%, this apparent discrepancy can be 
understood. 



3. Modeling the visibility data 

The most straight forward way of interpreting sparsely sam- 
pled interferometric data (visibilities) is by fitting the Fourier 
transform of an assumed brightness distribution of the object. 
Simple size estimations can be obtained by elementary models 
with only a few free parameters. As mentioned in the introduc- 
tion, the definition of a diameter is difficult because of its strong 
wavelength dependence, as well as of the intracycle and possi- 
ble long-term cycle-to-cycle variations (e.g. Hams' et al. 1 19951 1 
expected from models (e.g. Hofmann et al. I1998I I. On the other 
hand, the size and its wavelength-dependent shape can tell which 
layer of the atmosphere is actually observed and which chemi- 
cal and physical mechanism are responsible for this appearance. 
The low surface gravity results in an extended atmosphere and 
temperature structure and in the formation of molecular layers 
around late-type stars. Therefore, no sharp transition between 
star and circumstellar environment exists and limb-darkening ef- 
fects (center-to-limb variations) are very pronounced. 

All 75 visibility measurements of W Hya are plotted in Fig.|6] 
for one example wavelength bin (9.07 jjm) as function of spatial 
frequency (projected baseline divided by wavelength). There is a 
considerable spread in the data, but this should be expected since 
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Fig. 4. Same as Fig. [3] but plotted versus visual light phase. The MIDI observations are color-coded by position angle. The shaded 
and white regions (bins 1 to 4) ai^e used for intracycle studies in Sect. 14.21 



the measurements were obtained at different pulsation phases, 
pulsation cycles and position angles. In addition, some uncer- 
tainties in the reduction process remain. The left hand panel 
shows the data color-coded by visual light phase, while the right 
hand panel shows the data color-coded by position angle. In both 
plots systematically displaced distributions are noticeable. In the 
left hand panel, observations at visual maximum (dark points) 
have lower visibilities than those at pre-minimum (red points), 
if measurements with similar spatial frequencies are compared. 
However, this trend is very weak. In comparison, the displace- 
ment of the two distributions with different position angles is 
more significant in the right hand panel. For a discussion of this 
see Sect.m 

After the first zero, the visibilities in the second lobe at 
around 24 arcsec"' remain low. This indicates that a uniform 
disk (UD), with a constant brightness distribution up to the 
edge of the disk, cannot be applied to the data. A UD gives a 
higher second lobe in Fourier space. On the other hand, a simple 
Gaussian is not appropriate either, since a second lobe is clearly 
detected. Only a model that takes an extended atmosphere into 
account, which results in limb-darkening, can fit the data prop- 
erly. The simplest possibility with only a few free parameters and 
a low second lobe is a fully limb-darkened disk (FDD). To ac- 
count for a flux contribution of the extended silicate dust shell, 
the visibility function is not forced to be 1 at zero spatial fre- 
quency. Due to the lack of measurements at very low spatial 
frequencies, the dimension of the silicate dust shell cannot be 
constrained interferometrically. 



3. 1 . Fully limb-darkened disk (FDD) 

A fully limb-darkened disk model is fitted to the visibility mea- 
surements in order to get diameter estimates. The visibility and 



intensity functions of a circular FDD are given by 
3 ^/7T JiiiinOfBDr) 



V{r) = 



/(p) 



V2 (;r0FDDr)3/2 



and 



(1) 



(2) 







otherwise. 



respectively, where J2/2 is a Bessel function of the first kind of 
order 3/2, 0fdd the FDD diameter and e the flux contribution of 
the FDD. The flux contribution e can be less than 1 if a fully re- 
solved component (FRC), e.g. a surrounding silicate dust shell, 
adds flux at the given wavelength, since the total visibihty con- 
sists of both components: 



V, 



total 



kn^FDD+(l-e)n^, 



FRC I 



(3) 



The complex visibility of a fully resolved component, 'Vfrc, 
then does not contribute to the total visibility, but the flux 
fraction, /frc, of that component still modifies e via e = 
/fdd/(/fdd + /frc). The radius r in Fourier space is given by 
the spatial frequencies as r = Vm^ -i- u^. The spatial frequen- 
cies M and u are calculated from the projected baseline B, the 
position angle §, and the wavelength A via u = Bsin^/A and 
u - Bcos-d-fA. The radial coordinate p on the sky is defined 
by the angular coordinates a and /3 (angular separation rela- 
tive to the center of the object on the tangential sky plane) with 

p= yja^ +13^. 

The elliptical FDD, used in Sect. l4.1l to analyze asymmetries, 
is derived from the circular FDD by applying a rotation (with an 
orientation given by the position angle §) and a compression to 
one of the axes (which becomes the minor axis). The rotation can 
be obtained by changing the Fourier and sky reference frames 
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Fig. 5. Spectral energy distribution of W Hya. The included Photometry is not corrected for reddening effects and the black body 
curve is only included for guidance. The ISO LWS and SWS spectra from observation epochs 1996-02-07 and 1992-02-14 are used, 
respectively. The MIDI spectrum (green, where light green are the errors) does not show the silicate emission feature at around 
10 microns. The MIDI spectrum is shown in the inset in more detail. 



via 

M,j - u COS I? - y sin 
i>i> - M sin J? + t; cos -d- 

and 

a» - a cos J? - y6 sin 
P» - asmd- + p cos §, 



(4) 
(5) 

(6) 
(7) 



respectively. 

Applying the compression factor of 77 = ^FDD.min/^FDD.maj, 
with 0FDD,min the minor diameter and flpDD.maj the major diameter 
in Fourier space, yields the new variables 



rif,!! = ^^i^rf^^f, and (8) 
P.%r, = ^alli^^+pl (9) 
respectively. The visibility and intensity functions then become 

,maj 



V2 (Tr^FDD.maj'";*,,,)^^^ 



\ %DD.inai J 



and 

if 5< 1 
otherwise. 



(10) 



(11) 



respectively, with 



6 = 



fl2 

'^FDD.min 



'^FDD.maj 



(12) 



Best-model parameters are derived by performing the 
Levenberg-Marquardt least-squares minimization procedure 
programmed for the interactive data language IDL as MPFIT by 
C. B. MarkwardQ. To avoid the problem of running into local 
minima and to check for degeneracies, a global grid search was 
implemented. This is done in such a way that the least-squares 
minimization always operates only between two grid points in 
each parameter, while going through the whole parameter space. 
Formal 1 cr errors were computed by taking the difference be- 
tween the best fit value and the value given at the lowest non- 
reduced chi-square plus one, knowing that this is only an ap- 
proximation, since the chi-square values are not Gaussian dis- 
tributed. To get a more appropriate parameter error, the visibility 
errors are scaled to a value yielding a reduced chi-square of one. 
Besides this, the final parameter error is always a multiple of 
the shortest grid distance. This approach leads to an adequate 
description of the data by an FDD and to reasonable error esti- 
mates. 



3.2. Model fitting resuits 

The left hand panel of Fig. |6] shows the fit of a circular FDD 
and UD to all 75 visibility measurements for wavelength bin 6 



" http://cow.physics.wisc.edu/~craigm/idl/idl.html 
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Fig. 6. Fits of different models to the visibility measurements for wavelength bin 6 (9.07 //m). The 75 visibility data, taken over all 
pulsation phases and cycles, are plotted versus spatial frequency. Left: The measurements color-coded by visual light phase. The 
fit of a circular FDD gives for this wavelength a diameter of 79.2 mas and a relative flux contribution of 0.85. In addition, it can 
be clearly seen that a UD does not fit the data as well as an FDD. Right: The measurements color-coded by position angle. The 
elliptical FDD model fit results along the minor (76 mas) and major (85 mas) axis are included as solid lines. The fit to a given 
measurement is represented as a small open circle and is connected by a line with the corresponding observational point. The relative 
flux contribution is determined to be 0.86. 



(9.07 yum). Clearly, the second lobe is better fitted by an FDD. 
A comparison of the reduced chi-square values verified that a 
UD is not a good representation of the brightness distribution of 
W Hya. The best-fit parameters of the circular FDD are given in 
Table |2]for all 25 wavelength bins. Formal errors are computed 
as mentioned in the previous section, while the listed mean visi- 
bility errors are derived as described in Sect. 12.2] A residual anal- 
ysis shows that for each wavelength bin on average about 60% of 
the measurements are inside the 3 cr range. This low value stems 
from using rather low uncertainties for the visibilities and this fit 
not accounting for the scatter of the data due to the influences of 
the pulsation phase and the asymmetry. 

The left hand panel of Fig. Q shows the fit to all 25 wave- 
length bins. Clearly, the first zero shifts to lower spatial frequen- 
cies with increasing wavelength, i.e. the FDD diameter, 0fdd, 
increases with increasing wavelength. At the same time, the rel- 
ative flux contribution, e, of the FDD decreases. This behav- 
ior is shown more exactly in the two parameter plots in Fig. |8] 
From these panels it can be derived that ^fdd actually stays con- 
stant at a value of about (80 ± 1.2) mas between 8 and 10 /jm, 
while it gradually increases at wavelengths longer than 10 yum 
to reach (105 + 1.2) mas at 12 jum. The diameter increase in 
the longer wavelength regime corresponds to a relative increase, 
OFDD,\2tim/SFDD,\o^im, of (31 + 3)%. This apparent increase from 
80 mas to 105 mas is equivalent to an increase from 7.1 AU to 
9.5 AU at the distance of W Hya. The molecules and close dust 
species causing this shape are akeady indicated inside the plot, 
but will be discussed further in Sect. 13.41 In contrast, the relative 
flux decreases from (0.85 + 0.02) to (0.77 + 0.02), reflecting the 
increased flux contribution from the fully resolved colder sur- 
rounding silicate dust shell going to longer wavelengths. A fit to 
the unbinned grism data results in the same shapes and does not 
reveal any additional features. 

3.3. Wavelength dependence of the diameter 

Before interpreting the previously obtained FDD diameter be- 
havior in the next section, it is set in relation to measurements at 
other wavelengths. Figure |9] shows interferometric the diameter 



determinations reported by various authors from the visual to the 
mid-IR. They were obtained by fitting a Gaussian (green), a uni- 
form disk (UD, red) or a fully limb-darkened disk (FDD, blue). 
A conversion between the models is not performed, since the 
various diameter determinations depend on the spatial frequency 
coverageQ. Information on visual light phases and position an- 
gles (if applicable) of these observations can be found in Table[3] 
In addition to the FDD diameter, derived in this study from the 
fit to the full data set, the diameter trends at visual maximum 
(upper curve) and visual minimum (lower curve) are plotted as 
well (cf. Sect. 14:2b. 

The observed apparent diameter changes dramatically 
within the given wavelength range due to the strong 
wavelength-dependent opacities of the atmospheric constituents 
(cf. e.g. Baschek et al. 1991 and Scholz et al. 2001). In the op- 
tical, the measured diameters are sensitive to TiO bands. The 
largest variations are around the strongest bands at 712 nm and 
670 nm with apparent diameter enlargements of up to a factor 
of 2. It is notable that the diameters observed at around visual 
minimum by Ireland et al. (2004a) are systematically larger com- 
pared to the diameters obtained at around visual maximum by 
Tuthill et al. ( 1999). Ireland et al. (I2004al l conclude that signifi- 
cant dust production occurs near minimum, and suggest that the 
large increase in apparent diameter towards the blue is caused 
by light scattered by dust. Gaussian diameters between 0.7 and 
1 .0 fim were measured by Ireland et al. (i2004a) at two position 
angles (120° and 252°). Even though the diameters at both po- 
sition angles are slightly different, the authors conclude that this 
is not significant. 

In the near-IR, H2O and CO in different layers are predomi- 
nantly responsible for the wavelength dependence of the diame- 
ter. Most of the observations are conducted at J, H, and K bands 
(1.25, 1.65, and 2.2 fim, respectively). It can be inferred that also 
in the near-IR diameters vary by a factor of 2 and that there is 
a complex diameter dependence on pulsation cycle and pulsa- 



Empirical conversion factors are approximately: FDD » 1.15 UD 
x 1.68 FWHM (assuming a uniform spatial frequency coverage in the 
first lobe). 
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Table 2. Results of the circular and elliptical fully limb-darkened disk fits. 



circular FDD elliptical FDD 



^Vavelength 


Visibility 


fpDD 






Flux 6 


'7FDD,maj 


Axis Ratio t] 


PA& 


Flux e 


(jim) 


crror"^ 


(mas) 








(m&s) 






C) 








8.12 


0.007 


81.0 ± 


1 





0.86 


± 0.02 


92.5 ± 5.0 


0.87 


± 0.07 


-1 ± 


4 


0.89 ± 


0.03 


8.33 


0.007 


79.9 ± 


1 


2 


0.84 


± 0.01 


88.4 ± 5.0 


0.89 


± 0.08 


6 ± 


4 


0.87 ± 


0.03 


8.54 


0.006 


80.8 ± 


1 





0.83 


± 0.01 


89.3 ± 5.0 


0.89 


± 0.07 


-1 ± 


4 


0.84 ± 


0.03 


8.71 


0.006 


80.9 ± 


1 


6 


0.86 


± 0.02 


90.5 ± 5.0 


0.86 


± 0.07 


-1 ± 


4 


0.86 ± 


0.03 


8.87 


0.006 


79.6 ± 


1 





0.86 


± 0.01 


84.9 ± 5.0 


0.91 


± 0.08 


14 ± 


4 


0.86 ± 


0.03 


9.07 


0.006 


79.2 ± 


1 





0.85 


± 0.01 


85.3 ±5.0 


0.89 


± 0.08 


15 ± 


4 


0.86 ± 


0.03 


9.26 


0.006 


78.0 ± 


1 





0.83 


± 0.02 


85.1 ±5.0 


0.88 


± 0.08 


18 ± 


4 


0.85 ± 


0.03 


9.45 


0.006 


/ /.y ± 


1 
1 


2 


0.81 


± 0.02 


87.5 ± 5.0 


0.86 


± 0.08 


13 ± 


4 


0.83 ± 


0.03 


9.63 


0.006 


79.6 ± 


1 


2 


0.80 


± 0.02 


89.4 ± 5.0 


0.85 


± 0.07 


14 ± 


4 


0.81 ± 


0.03 


9.78 


0.006 


/9.8 ± 


1 
I 


2 


0.79 


± 0.02 


90.5 ± 5.0 


0.84 


± 0.07 


14 ± 


4 


0.81 ± 


0.03 


9.92 


0.006 


79.0 ± 


1 


2 


0.78 


± 0.01 


89.7 ± 5.0 


0.85 


± 0.07 


3 ± 


4 


0.78 ± 


0.03 


10.09 


0.005 


TOO 1 
10.0 ± 


1 
1 





0.77 


± 0.01 


89.7 ± 5.0 


0.85 


± 0.07 


3 ± 


4 


0.78 ± 


0.03 


10.26 


0.005 


80.4 ± 


1 





0.76 


± 0.01 


91.5 ±5.0 


0.85 


± 0.07 


5 ± 


4 


0.77 ± 


0.03 


10.42 


0.005 


84.5 ± 


1 


4 


0.77 


± 0.01 


93.7 ±5.0 


0.86 


± 0.07 


13 ± 


4 


0.77 ± 


0.03 


10.58 


0.005 


87.3 ± 


1 


4 


0.77 


±0.02 


98.1 ±5.0 


0.84 


± 0.07 


15 ± 


4 


0.78 ± 


0.03 


10.71 


0.005 


89.2 ± 


1 


4 


0.77 


±0.01 


99.6 ± 5.0 


0.85 


±0.07 


15 ± 


4 


0.78 ± 


0.03 


10.84 


0.005 


91.2 ± 


1 


2 


0.76 


±0.01 


102.2 ± 5.0 


0.86 


± 0.06 


11 ± 


4 


0.77 ± 


0.03 


10.99 


0.004 


94.8 ± 


1 


4 


0.76 


±0.01 


104.5 ± 5.0 


0.86 


±0.06 


15 ± 


4 


0.76 ± 


0.03 


11.14 


0.004 


96.3 ± 


1 


2 


0.76 


±0.01 


105.7 ± 5.0 


0.87 


±0.06 


14 ± 


4 


0.77 ± 


0.03 


11.29 


0.004 


97.8 ± 


1 


2 


0.77 


±0.01 


106.8 ± 5.0 


0.88 


±0.06 


15 ± 


4 


0.77 ± 


0.03 


11.43 


0.005 


98.9 ± 


1 


4 


0.76 


±0.01 


108.2 ± 5.0 


0.87 


±0.06 


15 ± 


4 


0.77 ± 


0.03 


11.55 


0.005 


100.5 ± 


1 


4 


0.77 


±0.01 


109.8 ± 5.0 


0.88 


±0.06 


14 ± 


4 


0.78 ± 


0.03 


11.67 


0.005 


101.6 ± 


1 


4 


0.77 


±0.01 


110.8 ±5.0 


0.88 


±0.06 


15 ± 


4 


0.78 ± 


0.03 


11.81 


0.004 


103.4 ± 


1 


4 


0.78 


± 0.02 


112.3 ±5.0 


0.88 


±0.06 


16 ± 


4 


0.79 ± 


0.03 


11.95 


0.004 


104.7 ± 


1 


2 


0.79 


± 0.01 


113.2 ±5.0 


0.88 


± 0.06 


19 ± 


4 


0.80 ± 


0.03 



This is the mean visibility error used for the corresponding wavelength bin. See Sect. l2.2l for the definition of this value. 




spatial frequency B / A (1/arcsec) position angie (deg) 



Fig. 7. Left: Same as left hand panel of Fig. |6] but for all 25 wavelength bins. Clearly the first zero shifts to lower spatial frequencies 
with increasing wavelength, i.e. the FDD diameter increases with increasing wavelength. Right: Visibility as function of position 
angle for wavelength bin 6 (9.07 /im). Measurements (crosses) are connected with their corresponding circular FDD model points 
(small open circles) showing that there is a systematic shift in the difference, Vmodei - ^'measurement, by going from 0° to 90°. For data 
points with low visibilities, corresponding to data in the increasing part of the second lobe, the systematic is inverse. See Sect.|4]for 
a detailed interpretation. 



tion phase (cf. e.g. Weiner et al. l2003l Woodruff et al. l2009b . 
The gap at around 2.7 yum in the curve of Woodruff et al. (2009) 
comes from telluric contamination. A reasonable estimate for 
the true photospheric diameter of an AGB star can be obtained 
from uniform disk measurements at K-Band. The 2.2 fim diam- 
eter is approximately 1.2 times the true photospheric diameter 
(e.g. Millan-Gabet et al. l2005T l. 



Wishnow et al. (1201 Oi l assume a UD diameter of the stel- 
lar component in the mid-IR (11.15 jum) of about 50 mas for 
their model. However, this value was reported as not very ac- 
curate. Since the FDD diameters obtained with MIDI describe 
a region whose exact location depends on the flux contribution 
of all constituents (continuum photosphere, atmospheric molec- 
ular layers, and nearby dust shells; cf. Sect. 13.4) as function of 
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Fig. 8. Left: The fitted circular fully limb-darkened disk diameter, ^fdd, as function of wavelength. The solid gray line shows the fit 
to all 75 measurements and the other lines the fits to individual pulsation phase bins. Right: The relative flux contribution, e, of the 
circular FDD obtained for the fit to the full data set. These values are fixed for the fit of the circular FDD to the individual pulsation 
phase bins. 



wavelength and pulsation phase, and not only the photosphere 
of the star, both diameters are not comparable, and the appar- 
ent discrepancy can be explained. The average of all diameters 
at 2.2 jum of the authors given in Table [3] is (42.8 + 3.5) mas. 
This gives a ratio, 6'FDD,iO/jm/6'uD,2.2Aim, of 1.8 + 0.2. The UD fit 
value should be used to compare similar models, and the ra- 
tio, 0uD,iO;/m/^^UD,2.2/im, bccomcs 1.6 + 0.2. From the mean K- 
band diameter the photospheric diameter, Ophoi, can be estimated 
to about 36 mas (42.8 mas divided by 1.2, cf. previous para- 
graph). This value is assumed for ^phot throughout this paper and 
gives a mid-IR to photospheric diameter ratio, 9\jd iO/jm/6'phot, of 
1.9 + 0.2. 

Maser observations of different molecules give additional di- 
ameter constraints. SiO and H2O masers probe inner regions 
where the molecular spheres are present and dust formation 
takes place, while OH masers trace wind regions farther out. 
Ring diameters for SiO masers were determined by Miyoshi et 
al. (fT994b . Cott on et a l. (l2004l 120081) . Reid & Menten (l2007l) . 
and Imai et al. (1201 Oi l. The average of their diameters is about 
(77 +11) mas, which gives a ratio, 0sio/%DD,iO/im, of 1.0 ±0.1 
(^sio/^UD.iO/im = 0.9 ±0.1). H2O maser ring diameters forW Hya 
are on the order of 300 mas (Reid & Menten ,1990 1. while OH 
maser ring diameters at 1665 and 1667 MHz were determined 
to about 700 and 1 130 mas, respectively (Szymczak et al. |1998l 
but see also Shintani et al. 2008 ). In addition, thermal line emis- 
sions of HCN (Mul ler et a l.[2008land Ziu rys et a l. [2009] l and CO 
(e.g. Olofsson et al. l2002l and Ziurys et al. |2009l l were repeatedly 
detected as well. 

3.4. Discussion and interpretation 

Studying the interplay of the extended pulsating atmosphere 
with molecular spheres directly above and the dust formation 
and wind acceleration zone farther out is crucial for understand- 
ing the physical and chemical processes working in AGB stars. 
With the MIDI interferometer, the photosphere and molecular 
layers of W Hya are probed, sampling also the region where the 
first seeds for dust formation originate. 

The spectrum in Fig. |5] shows that molecules, such as CO, 
H2O, and SiO, are present in the upper atmosphere. In the N- 
band between 8 and 1 3 fim, strong pure -rotation lines of H2O are 



expected. In addition, SiO exhibits fundamental bands between 
8 and 10 yum (Decin 12000b . Modeling by Justtanont et al. ( 120041 ) 
also reveals at least three different dust species, namely amor- 
phous silicate, AI2O3, and MgFeO, which are likely located at 
different distances from the star. In particular, amorphous AI2O3 
provides significant opacity for wavelengths longwards of 1 jum 
(cf. e.g. Koike et al. 1995, Begemann et al. 1997, Posch et 
al. 1999 Egan & Sloan Woitk e et a l. 11006 Ireland & 

Scholz 2006 and Robinson & Maldoni [2010l) . and it can survive 
at high temperatures. 

Quasi-static, warm, and dense molecular layers close to the 
star, at typically 2-3 photospheric radii (Rphot), are detected for 
0-rich (e.g. Mennesson et al . 120021 Perrin et al. 120041 Ireland 
et al. I2004dl Woodruff' et al. EOOTand Fedel e et al . '2005 ) and 
C-rich (e.g. Hron et al. . 19981 and Ohnaka et a l 12007^ AGB stars, 
as weU as for RSG stars (e.g. Perrin et al. 12007) . These lay- 
ers were introduced earlier to explain spectroscopic observations 
(e.g. Hinkle & Barnes 1979; Tsuji et al. 1 19971 and Yamamura et 
al. lT999] l. 

In particular, the O-rich Miras RR Sco (Ohnaka et al. 2005) 
and S Ori (Wittkowski et al. 2007) show a diameter behavior 
very similar to that of W Hya throughout the N-band. Comparing 
the shape from the left hand panel of Fig. |8] with Fig. le in 
Ohnaka et al. (2005 ) and Figs. 2 - 5d in Wittkowski et al. (l2007l l 
leads to the idea that, in all three stars, the same constituents and 
mechanisms are responsible for this appearance. 

The diameter ratio 0i2/jm/SiO/jm is approximately 1.3 for 
RR Sco, for S Ori between 1.3 and 1.5, and for W Hya ap- 
proximately 1.3 (Sect. [JST l. Even the increase between the K- 
band diameter and the diameter at 10 /im is similar. The ratio 
6'uD,iO/jm/^^UD,2.2//m is for RR Sco approximately 1.8, for S Ori be- 
tween 1.6 and 2.2 and for W Hya approximately 1.6 (Sect. l33T l. 

Ohnaka et al. ( 120051 ) modeled RR Sco by adding a molecular 
layer of H2O and SiO, and a dust layer, consisting of aluminum 
oxide and silicates, to a static star with a fixed temperature. This 
composition reproduced the N-band diameter behavior very rea- 
sonably, with an inner radius for the dust shell of 7.5 .Rphot, while 
the molecular layers were located at a radius of 2.3 Rphot- The 
large dust shell radius might be due to the mix of 80% aluminum 
oxide (AI2O3) and 20% silicates and because the same density 
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Fig. 9. Diameter measurements of W Hya over a wide wavelength range with different models reported by various authors. In 
addition to the N-band full data set diameter (thick blue line; all 75 observations), the diameters at maximum light (upper thin blue 
line) and minimum light (lower thin blue line) from this work are shown as well (FDD » 1.15 UD x 1.68 Gaussian). See Table [3] 
for more details. 



Table 3. Details to the diameter measurements presented in Fig. |9] 





Author 


Model 


Phase" 


Comments and position angles (PA) 


A 


Lattanzi et al. f 1997)1 


UD 


0.64 


two perpendicular axes (rj » 0.86, PA a 143°, 583 nm) 


B 


Ireland et al. ( 20Q4all 


Gaussian 


0.44 


upper curve for PA 120° and lower curve for PA 252° (0.7 - 1.0 iJ.m) 


C 


Haniff etal. (1995) 


UD 


0.04 


uniform disk diameter at 700 and 710 nm 


D 


Tuthill etal. (199^1 


Gaussian 


0.04 


elliptical Gaussian at 700 and 710 nm (77 « 0.94, PA a 93°) 


E 


Monnier et al. (12004b 


UD 


0.50 


UD model was reported as a bad fit to the data (K band) 


F 


Woodruff et al. 12009J 


UD 


0.58- 1.53 


curves for phase 0.58 (middle), 0.79 (lower) and 1.53 (upper) 


G 


Woodmff et al. d2008k 


UD 


0.50- 1.00 


mean of multiple measurements in given phase range (I, H, K, & L band) 


H 


Millan-Gabet et al. (12005b 


UD 


0.60 


UD diameter for H and K bands 


I: 


Wishnow et al. POIOJ 


UD 


0.1 - 1.1 


assumed photospheric UD diameter at 11.15 /im, reported as not very accurate 


J: 


this work 


FDD 


0.00 & 0.50 


curves for full data set (middle), max (upper), and min (lower) light 



visual light phase (phases refer to the phases determined by the original authors) 



profile and inner condensation radius were used for both dust 
species. 

Wittkowski et al. (I2007I I modeled S Ori with the dynamic 
atmospheric models of Ireland et al. (I2004bl I2004cb . where the 
molecular layers are naturally included and only the dust shell 
has to be added ad hoc. This assumes that the stellar atmosphere 
and dust shell are spatially separated. They allowed different 
density profiles and condensation radii for the aluminum ox- 
ide and silicate dust. The model without silicates fitted the data 
better. The inner boundary of the AI2O3 dust shell varied be- 
tween 1.8 and 2.4 photospheric radii. At visual minimum, the 
dust shell was more compact with a larger optical depth, and it 
had a smaller radius, while at visual maximum the opposite was 
the case. 

W Hya can be explained by a similar model. The overall 
larger diameter in the mid-IR is caused by a warm molecular 



layer of H2O, and the gradual increase longward of 10 //m arises 
most likely from the presence of AI2O3 dust close to this layer. 
A nearby SiO molecular shell might be of some relevance for 
the diameter enlargement between 8 and 10 /im as well. The 
formation of AI2O3 dust at these short distances from the stellar 
surface would be consistent with the empirical results by Lorenz- 
Martins & Pompeia ( I2000I I. as well as with recent theoretical 
calculations by Ireland & Scholz (2006) and Woitke (2006). 

Dust can only exist close to the star if the temperature is be- 
low the condensation temperature of the dust species. A rough 
estimate of the equivalent blackbody temperature from the to- 
tal MIDI flux, the flux fraction e, and the diameter 0fdd gives 
a value of about (800 ± 100) K at the location of the FDD ra- 
dius. This temperature is low enough that aluminum oxide dust 
can condensate (cf. e.g. Woitke II999I I. Since AI2O3 has only 
a moderate abundance and a low absorption efficiency at op- 
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tical and near-IR wavelengths, i.e. near the flux maximum at 
around 1 jum, it cannot be responsible for initiating the mass loss 
(cf. e.g. Woitke 2006). AI2O3 can exist close to the star without 
inducing mass loss. 

Scattering off large Fe-free silicate grains may solve this 
problem, as these micron-sized silicates do not have a con- 
siderable absorption cross-section, but a high radiative scat- 
tering cross-section around the flux maximum at 1 fim 
(cf. Hofner l2008l l. However, these grains are not detected with 
MIDI. If a large amount of Fe-free silicates existed close to the 
star, an emission feature in the MIDI spectrum at around 10 fim 
would be present (due to strong vibrational resonances), and it 
would also modify the apparent FDD diameter around this wave- 
length. Therefore, the wind acceleration mechanism proposed 
by Hofner (2008) cannot be directly verified. However, since 
the mass-loss rate of W Hya is very low (cf. e.g. Nowotny et 
al. 120 101) . this behavior is not surprising. 

This applies not only to W Hya, but also to RR Sco and S Ori. 
All three O-rich AGB stars have moderate to low mass-loss rates. 
The similar behavior of the apparent diameter throughout the N- 
band suggests that the appropriate mechanism for a moderate-to- 
low mass-loss rate is similar, since the conditions appear similar 
in the transition region (between the photosphere and the sili- 
cate dust shell) where the wind should be initiated, even if the 
mechanism cannot be clarified here. 

In general, another possibility of forming a wind in O- 
rich AGB stars is to accelerate small amounts of carbon grains 
(Hofner & Anderson I2007t . Unfortunately, these carbon grains, 
formed in nonequilibrium environments, do not show any spec- 
tral features in the mid-infrared and would probably produce IR 
colors that are not consistent with the observations (Hofner, pri- 
vate communication). In this context it might be interesting to 
know whether scattering on AI2O3 grains is important, even if 
they are probably not abundant enough. Also the role of large 
amounts of water vapor in molecular shells and the radiation 
pressure on water molecules may need more detailed calcula- 
tions. 

As described in the previous section, W Hya exhibits SiO 
maser emission at a distance of about 0.9 times the equivalent 
UD diameter. The SiO maser spots are therefore co-located with 
the close aluminum oxide dust shell. This is very similar to what 
was found for S Ori (Wittkowski et al. I2007K and supports an 
explanation by a similar model. 

From the analysis of the spectrum and the fact that the flux 
contribution of the star/layer is below one, it was concluded that 
a circumstellar silicate dust shell exists around W Hya with a 
lower limit to the condensation radius of about 28 times the 
photospheric radius. At these large radii, aluminum oxide dust 
grains are no longer detectable, since silicon and magnesium are 
much more abundant than aluminum, and they dominate the dust 
opacities. The low flux contribution of the silicate dust shell is 
consistent with a low mass-loss rate of W Hya. 

With the previous results, the shape of the visibility curves in 
the left hand panel of Fig. [T]can be understood qualitatively. The 
visibility increase between 8 and 10 fj.m corresponds to the par- 
tially resolved stellar disk (including the close molecular layers) 
at an almost constant FDD diameter, while in the region between 
10 and 12 yum, the flux contribution of the spatially resolved 
AI2O3 dust shell becomes more relevant. This leads longward 
of 10 ;um to the nearly flat visibility curve and an increased FDD 
diameter (cf. Ohnaka et al. 2005 and Wittkowski et al. 2007). In 
addition, extinction of the outer silicate dust shell becomes im- 
portant. The comparison with RR Sco and S Ori suggests that 



the partially resolved molecular layers are optically thick and 
that the nearby AI2O3 dust shell is optically thin. 

It is clear that quantitative modeling is needed to support the 
above findings, in particular, if the derived constituents of the 
close molecular and dust layers (H2O, SiO?, AI2O3) can really 
provide sufficient opacities to explain the observed diameter de- 
pendence on wavelength, in particular if AI2O3 could cause the 
apparent diameter increase beyond 10 jim. To quantify the re- 
sults, it will be necessary to apply dynamic atmospheric models 
of e.g. Ireland & Scholz (2006 ). Even if AI2O3 dust has not been 
consistently included in these models so far, and has often to 
be added ad hocQ This will give more detailed insight into the 
physical processes at work there. 

4. Asymmetry, intracycle variations, and 
cycle-to-cycle variations 

Since observations at similar pulsation phases were mostly con- 
ducted at similar position angles, the effects of different diam- 
eters, owing to asymmetry and pulsation, are unfortunately not 
easy to disentangle. At visual minimum, the position angles (PA) 
cluster at around 90° (green tick marks in Fig.|4}. At these posi- 
tion angles, the visibilities are higher than expected for a circular 
FDD, i.e., a smaller diameter is observed. At visual maximum, 
the PA is around 50° (red tick marks in Fig. |4) with the result of 
getting lower visibilities as expected for a circular FDD, i.e., a 
larger diameter. This trend is more clearly shown in the right 
hand panel of Fig. |7] where the visibilities of the observations 
and the visibilities of a circular FDD are compared as a func- 
tion of position angle. Obviously, the difference between both 
quantities changes its sign by going from 0° to 90°. 

The above behavior is notable if all 75 observations over all 
three pulsation cycles are included in the investigation. A care- 
ful analysis, as described below, reveals both effects in W Hya. 
The result will be that the diameter variation due to an elliptical 
asymmetry and due to the pulsation effect have the same order 
of magnitude. 

The pulsation cycle is divided into four bins as shown in 
Fig-El Bins 1, 2, and 3 consist of 23, 42, and 10 observations 
with phase ranges 0.875-0. 125, 0.125-0.375, and 0.375-0.625, 
respectively. There are no observations for bin 4 with a phase 
range of 0.625 - 0.875. Bin 1 contains observations at visual 
maximum and bin 3 observations at visual minimum. 

4.1. Elliptical asymmetry 

Studying asymmetric features requires a good angular uv- 
coverage at each projected baseline length obtained within a rea- 
sonable time. Since the dependence on the visual light phase in- 
troduces akeady additional complexity into the model, only a 
simple elliptical model can be applied in order to get an indica- 
tion for a departure from symmetry. 

For these reasons, an elliptical FDD is fitted to subsets of the 
full data set. Each subset contains data with similar spatial fre- 
quencies and a narrow visual light phase range at a certain pul- 
sation cycle. For all these configurations the resulting fits show 
similar nonuniformity with comparable position angles and axis 
ratios. The fit was repeated for the full data set to obtain an over- 
all departure from symmetry and to have a more reliable uncer- 
tainty estimation, since the subsets have a high scatter and can 

Dynamic atmospheric models, including both the molecular layers 
and close dust formation, are in development (Scholz and Ireland, pri- 
vate communication). 
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Fig. 10. Asymmetric model of W Hya. Shown are the elliptical FDD fits for wavelength bin 6 (9.07 //m). Left: The model visibility 
(yellow - high visibility, blue = low visibility) plotted together with the data points in the uv-plane. The major axis in Fourier space 
has a position angle of 105°. Right: The corresponding image with a position angle of the major axis of 15° on the sky (linear 
intensity scale in arbitrary units). North is up and east is left. The size is 150 mas x 150 mas. 



contain discordant values due to fewer data points. The similarity 
of the results for the subsets shows that the ellipticity is stable to 
first order over a few years, over a certain stellar extension, and 
not dependent on the pulsation phase. Therefore, the parameters 
obtained for the full data set are given in Table |2] 

Listed are the major diameters, flpDCmaj, the axis ratios, rj, 
the position angles, i?, and the relative fluxes, e. The mean i? 
and T] over the full wavelength range are (11.2 ± 6.2)° and 
(0.87 + 0.07), respectively. Compared to the symmetric FDD, 
the fluxes are identical within the errors as expected, while the 
reduced chi square estimates gave lower values, thus indicating 
a slightly better model. The right hand panel of Fig.|6]shows the 
fit of the elliptical FDD for wavelength bin 6 (9.07 /^m). The 
visibility model curves for the minor and major axes are drawn 
together with the data points, showing that they fit the measure- 
ments at the corresponding angles. 

However, the parameter errors are most probably underesti- 
mated because of the high number of free parameters in com- 
bination with not negligible uncertainties in the data. In par- 
ticular, the PA estimation is not very well constrained, and the 
error might be more on the order of 15° to 20°, if compared 
with all fitted subsets. Figure [10] shows the asymmetric appear- 
ance of W Hya. The left hand panel shows the fit of the ellip- 
tical FDD for a representative wavelength, and the right hand 
panel displays the corresponding appearance on the sky for this 
model and mid-lR wavelength. The investigation of the differen- 
tial phases obtained with MIDI are unfortunately not conclusive, 
since they are smaller than the assumed errors, and no additional 
constrains can be inferred from them. 

Position angle estimations in the literature are very contra- 
dictory and are summarized in Tableland illustrated in Fig. [TT| 
Lattanzi et al. (.1997^ 1 and TuthiU et al. ( 1999 ) have found in the 
visual a PA of 143° and (94 ± 33)°, respectively, while no asym- 
metries could be detected in the near IR within the measure- 
ment uncertainties (Ireland et al. 2004a and Monnier et al. 2004). 
From the photospheric radio continuum, Reid & Menten (I2007t 



Since observations were obtained only on two perpendicular base- 
lines, the true major axis cannot be recovered. 



determine a PA of (83 ± 18)°, whereas Szymczak et al. (1 19981 1 
find for the projected angle of the magnetic field on the plane of 
the sky, obtained from OH maser observations, a value of about 
-20° (the PA of the underlying OH maser distribution is 70°), 
and a velocity discontinuity along the N^S axis. The 215.2 GHz 
SO line, observed by Vlemmings et al. ( 1201 I V also shows an ve- 
locity gradient along an N^S axis with an PA of the structure 
of (3 + 10)°. MuUer et al. (2008) detect in the HCN velocity 
an NW— >SE gradient. Interestingly, the outer dust envelope, ob- 
served by Marengo et al. (I2000I I at 18 yum, shows a very similar 
N-S elongation as determined in this work. Asymmetries in the 
SiO and H2O maser spot distributions were, e.g., reported by 
Reid & Menten ( 1 1990120071) and Imai et al. (l20T0b . 

From Fig. [TT| it can be inferred that the reported velocity 
gradients and the magnetic field alignment are approximately 
in a direction perpendicular to an E-W axis where most of the 
maser spots are located. If the velocity gradients originate from 
a general rotation of the circumstellar environment and the maser 
spots trace the highest molecule concentrations, then this might 
suggest that molecular gas is primarily ejected along this E-W 
polar axis. This would be also consistent with the PA of the 
photospheric radio continuum. The dust emission traced by the 
observations at 18 yum and the N-band are approximately elon- 
gated in a more perpendicular direction to this. This might sug- 
gest that the dust is more concentrated in an equatorial disk (or 
ring) located along an N-S axis. Therefore, one could conclude 
that the molecular gas and the dust are spatially separated, i.e., 
they have preferred ejection directions, and the different obser- 
vational wavelengths probe this. 

However, this is still very uncertain since not all PA determi- 
nations fit into this picture or deviate quite a lot from this, e.g., 
the PAs at optical wavelengths and the 215 GHz PA. Therefore, a 
strict connection between the outer structures and the close stel- 
lar regions also cannot really be established. One has to keep in 
mind that the concentration of maser spots along the E-W polar 
axis may only reflect that no velocity gradients exists along this 
axis and masers are traceable because of the inexistent Doppler 
shifts. 
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Fig. 11. Asymmetry and maser measurements of W Hya reported by various authors. See Table|4]for references and text for expla- 
nation. 

Table 4. Position angle (PA) determinations by various authors (see Fig.[TT]for an illustration). 



Author 


Ref" 


Obs Date 


PAC>) 


Wavelength 


Region'' 


Lattanzi etal. (11997b 


[8] 


1995 Dec 17 


143" 


583 nm 


inner 


Tuthill et al. (119991 


[7] 


1993 Jun 


94 ±33 


700 and 710 nm 


inner 


Ireland et al. (I2004at 




2001 Feb 8/9 


_d 


680 - 940 nm 


inner 


Monnier et al. (I2004J 




2000 Jan/Feb 


_d 


2.25 iim (K band) 


inner 


this work 


[9] 


2007 - 2009 


11 ± 20 


8 - 12 //m (N band) 


inner 


Reid & Menten j2007|l 


[6] 


2000 Oct/Nov 


83 ± 18 


43 GHz (radio continuum) 


inner 




[5] 


2000 Oct/Nov 


contours 


43 GHz (SiO maser) 


inner 


Imai et al. ( 2010) 


[10] 


2009 Feb 


spots 


43 GHz (SiO maser) 


inner 


Reid & Menten i'l990ll 


[4] 


1990 Feb 


contours 


22 GHz (H2O maser) 


intermediate 


Marengo et al. (2000 ) 


e 


1999 Jun 


N-S 


1 8 yum ( 1 . 1 8"x 1 .45" FWHM) 


outer 


Szymczak et al. (|1998|l 


[1] 


1996 Jan 15 


N-^S 


1667 MHz (OH maser), velocity gradient 


outer 




[2] 


1996 Jan 15 


-20 ± 20^ 


1667 MHz (OH maser), magnetic field 


outer 


MuUeretal. (12008 k 


[3] 


2008 Apr 13/15 


NW-»SE 


266 GHz (thermal HCN), velocity gradient 


outer 


Vlemminas et al. (12011b 


e 


2008 Jul 20 


3+10 


215 GHz (SO 55 - 44) 


outer 



" numbering used in Fig.[TT] ^ probed regions: inner = 50 - 100 mas, outer a 1000 mas; the true major axis cannot be recovered since 
observations are at two perpendicular baselines; no departures from symmetry could be detected within the measurement uncertainties; not 
shown in Fig.[TT]since the whole structure does not fit into the illustration; ' the position angle of the major axis of the underlying OH maser 
distribution is 70° 



Other explanations are therefore still possible, in particu- 
lar if time-dependent effects are important. The wide range of 
measured PAs at close photospheric radii might be related to a 
nonradial pulsation or due to clumpy dust formation. A tempo- 
rally varying non-symmetric brightness distribution could also 
be caused by stellar spots or large convection zones. A defor- 
mation as a result of a close companion can most probably be 
excluded, since no evidence for such a scenario has been ob- 
served so far. However, from the data collected here, no signif- 
icant temporal trends could be derived. In particular, the MIDI 
observations did not show any time variations of the asymmetry 
over about three years. This calls for more high-resolution obser- 
vations and more detailed modeling. A connection to the highly 



elongated structures seen in planetary nebulae cannot be drawn, 
since W Hya is believed to be still in an early AGB phase. 

4.2. Intracycle variations 

To estimate a light phase dependent angular diameter in the pres- 
ence of a position angle dependence, the projected baseline val- 
ues of the input data are transformed in a way such that the el- 
liptical model transforms into a circular model. This means that 
the projected baselines are artificially shortened (around the ma- 
jor axis) or lengthened (around the minor axis), i.e. as function of 
position angle, so that an elliptical model fit to the data shows no 
departure from circular symmetry in the new coordinate system. 
The parameters derived from the elliptical FFD fit are used for 
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this shearing. This transformation keeps the fitted absolute di- 
ameter values from being meaningful, so only differences with 
respect to a fit of a circular FDD to the full sheared data set are 
given. 

The fit of a circular FDD to the sheared data of the inter- 
mediate phase bin (bin 2) gives on average a (1.2 + 2.3) mas 
smaller diameter. The difference between maximum phase, with 
an on average (2.2 ± 2.3) mas larger diameter, and minimum 
phase, with an on average (3.2 + 3.7) mas smaller diameter, is 
(5.4 + 1.8) mas (cf. left hand panel of Table|5]l- This corresponds 
to a percentage change of (6 + 2)%. The individual wavelength- 
dependent trends can be seen in the left hand panel of Figure |8] 
As can be inferred from the plot, the shapes for all individual 
phase bins are very similar, indicating again that the different 
constituents probed at different wavelengths behave similarly. 
The relative flux was fixed to the value obtained for the full data 
set (shown in the right hand panel of Fig.O. 

The diameter variation is relatively small compared to the 
findings by Wittkowski et al. (I2007I I for the Mira variable S Ori. 
A phase-dependent diameter of S Ori cannot be derived directly 
from their measurements, since observations at visual maxi- 
mum and visual minimum were executed at different projected 
baselines. However, from the modeling results shown in their 
Fig. 12, it can be roughly estimated that the photospheric diam- 
eter changes by approximately (15 ± 5)% and show a similar 
behavior as function of visual phase. The smaller variation for 
W Hya can be expected since it has a smaller visual amplitude 
and is classified as an SRa/Mira variable, but it could also be 
related to the large phase-binning used in this study. 

The observed smaller angular diameter at visual minimum 
and the larger angular diameter at visual maximum can be 
explained by the phase-dependent presence of water vapor 
(Matsuura et al. 120021) and aluminum oxide dust (Ireland & 
Scholz 120061 and Witflcowski et al. 120071) . The star is hotter at 
visual maximum and H2O and AI2O3 dust can only exist farther 
out. At visual minimum, dust forms closer to the stellar surface 
and in larger amounts (otherwise it would be not detectable with 
MIDI close to the star), and suggests that the mass-loss rate is 
higher in this phase or shortly after this phase. 

While the observations made in this work show that possible 
AI2O3 dust exists closer to the star at visual minimum and farther 
out at visual maximum, it cannot be determined conclusively 
how this relates to the acceleration of the wind and the low mass- 
loss rate. AI2O3 can exist in the upper atmosphere without mass 
loss, and potentially relevant micrometer-sized Fe-free silicates 
were not detected with MIDI due to the low mass-loss rate and 
a therefore low abundance. Therefore, it is still possible that the 
wind formation and mass-loss mechanism is to a certain amount 
decoupled from the pulsation, even if it can be speculated that 
the relevant constituents for the wind acceleration are formed 
along with AI2O3. 

Simultaneous observations, tracing the relevant constituents 
and providing kinematic information, are necessary (e.g., high- 
resolution radio and mm observations of molecules/masers), and 
a comparison with dynamic atmospheric models can only be of 
limited success as long as the dust formation is not included. 

4.3. Cycle-to-cycle variations 

W Hya was observed over three consecutive pulsation cycles. 
For each cycle, data between phases 0.2 and 0.4 are used 
to search for cycle-to-cycle variations. These three ranges are 
shaded in Fig. |3] The ranges contain 8, 8, and 1 1 observations, 
respectively, and have a sufficient spatial frequency coverage. A 



circular FDD is fitted to these observations (which are sheared 
as described in Sect. 14.2b . while the relative flux is fixed to the 
value obtained for the full data set. This compensates for fewer 
measurements and avoids the problem of being too dependent 
on fewer points at low spatial frequencies (i.e. < 10 arcsec"'). 

Cycles 1 and 2 give similar diameters, with values on average 
of about (3.0 ± 2.4) mas below the diameter obtained for the full 
data set (cf. right hand panel of Table|5]l. The diameter for cycle 3 
is marginally higher, being on average (1.5 ± 3.0) mas above the 
full data set value. This means that the maximal variation due to 
non-repeatability of the pulsation is on the order of (5 ± 4)%, 
and is thus lower than intracycle variations. However, regard- 
ing the uncertainties, this is only marginally significant. The 
wavelength-dependent shape of the diameter does not change 
notably between cycles. 

This gives an order of magnitude estimate of how large 
pulsation instabilities in SR variables are, but also reflects 
the chaotic atmospheric behavior. These variations are smaller 
than might be expected. From models by e.g. Hofner & 
Dorfi ( 1997 ), Hofmann et al. 1998 Ireland et al. (2004c), Ireland 
& Scholz (2006), and Nowotny et al. (2010) for Mira variables, 
it can be derived that the location of mass shells and dust con- 
densation radii can vary by up to 20% with a characteristic mean 
of 10%. However, it should be kept in mind that the movement 
of mass shells and dust condensation radii cannot be directly 
compared to changes in the atmospheric radius-density structure 
that is traced by interferometric observations. On the other hand, 
these marginal size changes are not surprising if compared with 
the relatively small intracycle variations and the fact that W Hya 
is an SR variable. Finally, this verifies that the folding of con- 
secutive pulsation cycles in the previous diameter and intracycle 
analysis is an acceptable assumption. 



5. Summary 

W Hya was monitored over about three years in the thermal IR 
(8 - 12 fj.m). These are the first high-resolution interferometric 
N-band observations of W Hya with MIDI. A photometric study 
reveals a clear phase dependency of the N-band flux with a flux 
variation on the order of 20% between maximum and minimum 
light. The mid-lR maximum occurs after the visual maximum at 
visual phase 0. 15 + 0.05. 

The visibility data can be best fitted with a fully limb- 
darkened disk, which accounts to some extent for surround- 
ing atmospheric layers. The resulting FDD diameter of W Hya 
is almost constant between 8 and 10 fim at a value of about 
(80 + 1.2) mas (7.8 AU), while it gradually increases at wave- 
lengths longer than 10 //m to reach (105 ± 1.2) mas (10.3 AU) at 
12 fj.m. In contrast, the relative flux decreases from (0.85 ± 0.02) 
to (0.77 + 0.02), reflecting the increased flux contribution from a 
fully resolved surrounding silicate dust shell. From field-of-view 
effects, it could be derived that the silicate dust shell has an inner 
radius larger than 28 photospheric radii (>50 AU, >0.5 arcsec). 

The measured apparent mid-IR diameter at 10 /im is 
about 1.6 + 0.2 times larger than the near K-band diameter 
(^^UD,iO/jm/6'uD,2.2/im), i-C- about 1.9 times the photospheric diame- 
ter. This is very similar to findings by Weiner et al. (120031 . Perrin 
et al. (120041 ) and Tatebe et al. (120061 ) for similar stars. In particu- 
lar, the diameter behavior throughout the N-band is comparable 
with observations of the O-rich Mira stars RR Sco and S Ori by 
Ohnaka et al. (2005 ) and Wittkowski et al. (2007 ), respectively. 

W Hya is therefore described by an analogous model. The 
constant diameter part results from a partially resolved stellar 
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Table 5. Intracycle (left) and cycle-to-cycle (right) variations. 



Phase 


Phase 


No of 


Absolute" FDD 


Relative" FDD 


Cycle'' 


No of 


Absolute" FDD 


Relative" FDD 


Bin" 


Range 


Obs. 


Difference (mas) 


Diameter 


Obs. 


Dilference (mas) 


Diameter 


1 


0.875-0.125 


23 


+(2.2 ± 2.3) 


1.02 ± 0.02 


1 


8 


-(2.9 ± 2.4) 


0.97 ± 0.02 


2 


0.125-0.375 


42 


-(1.2 ±2.3) 


0.99 ± 0.02 


2 


8 


-(3.1+2.4) 


0.96 ± 0.02 


3 


0.375-0.625 


10 


-(3.2 ± 3.7) 


0.96 ± 0.04 


3 


11 


+(1.5 ±3.0) 


1.02 ± 0.03 


4 


0.625-0.875 


















" cf. Fig.|4j cf. Fig.[3j " in respect to the full data set value 



disk, including the close molecular layer of H2O, while the in- 
crease beyond 10 /urn can be most likely attributed to the con- 
tribution of a spatially resolved nearby amorphous AI2O3 dust 
shell. Probably owing to the low mass-loss rate, close Fe-free 
silicate dust, as proposed by Hofner 2008 could not be detected 
with MIDI. 

Since observations at similar pulsation phases were mostly 
conducted at similar position angles, the effects of different di- 
ameters, owing to elliptical asymmetry and pulsation, are unfor- 
tunately not easy to disentangle. By only using certain pulsation 
phases and projected baselines an asymmetric character of the 
extended structure could be confirmed. An elliptical FDD, with 
a position angle of (11 + 20)° and an axis ratio of (0.87 + 0.07), 
fits the data. The asymmetry might be explained by an enhanced 
dust concentration along an N-S axis. 

To estimate the dependence of the angular diameter as 
a function of visual light phase, the input data are sheared. 
The observed angular diameter is smaller at visual minimum 
and larger at visual maximum with a periodic change of 
(5.4 + 1.8) mas between maximum and minimum, coiTespond- 
ing to about (6 + 2)%. This is much less than reported for 
the high amplitude pulsating Mira star S Ori (Wittkowski et 
al.p007). The smaller observed angular diameter at visual min- 
imum can be explained by the phase-dependent presence of wa- 
ter vapor and likely aluminum oxide dust and their temperature 
sensitivity. Since this variation only traces the location of con- 
stituents that are probably not relevant for the wind formation, no 
firm conclusions can be drawn concerning the mass-loss mecha- 
nism. One can only speculate that more dust forms at visual min- 
imum and that the mechanism for a moderate-to-low mass-loss 
rate is similar for 0-iich SR and Mira variables. The detected 
cycle-to-cycle variations are smaller than intracycle variations 
and on the order of (5 ±4)%. 

The observation of nearby molecular layers and a nearby 
dust shell, like in other AGB stars, confirms the emerging stan- 
dard picture and supports the need for self-consistent dynamic 
atmospheric models with consistently included dust formation 
close to the star. Primarily, the close AI2O3 dust shell most 
likely detected in our observations, well below the distance at 
which the silicate dust shell is traced, has now been revealed in 
a few objects. It has also been shown that a good phase and uv- 
coverage over the whole pulsation cycle can be crucial, while 
interferometric observations in the N-band are an irreplaceable 
tool for resolving close stellar structures and for searching for 
atmospheric constituents. Future work will concentrate on im- 
proving the data reduction process further and on applying dy- 
namical atmospheric models to the data. 
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